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ABSTRACT 

From a survey of the 3400 A region in the earliest 0-type spectra, we have fonnd that 
two of the four 02 giants observed in the Large MageUanic Cloud have the O IV lines 
there stronger than the N IV, while the other two have the opposite. A Small Magellanic 
Cloud counterpart also has the N IV stronger than O IV. Inspection of the blue spectra 
of these stars shows that the former pair have weaker N lines in all ionization states (III, 
IV, V) present as well as lines of C IV A4658, while the latter three have stronger N lines 
and greater He/H. Space ultraviolet observations of two of the N-strong stars show N V 
wind profiles substantially stronger than those of C IV, while in the N-weak stars the C IV 
features are equal to or stronger than the N V. The N-strong stars are now reclassified as 
0N2 III(f*), newly defining that category. These characteristics strongly suggest a larger 
fraction of processed material in the atmospheres of the 0N2 stars, which we confirm by 
modeling the optical spectra. In the context of current models, it is in turn implied that 
the 0N2 stars are in a more advanced evolutionary state than the others, and/or that they 
had higher initial rotational velocities. The recent formulation of the effects of rotation on 
massive stellar evolution introduces an additional fundamental parameter, which the CNO 
abundances are in principle able to constrain. We present some illustrative comparisons with 
current Geneva evolutionary models for rotating massive stars. It is possible that these very 
hot, nitrogen-rich objects are products of homogeneous evolution. Our results will provide 
motivation for further physical modeling of the atmospheres and evolutionary histories of 
the most massive hot stars. 

Subject headings: Magellanic Clouds — stars: abundances — stars: early-type — stars: 
evolution — stars: fundamental parameters 



1. INTRODUCTION 



Anomalies in the relative intensities of CNO lines 
in OB spectra have been described for some time 
(Walborn 1976, 1988, 2003; and references therein). 
Anticorrelations of N vs. C, O and correlations with 
He/H in the expected sense have encouraged inter- 
pretations in terms of mixing of CNO-cycled ma- 
terial into the atmospheres and winds of massive 
stars. An important point has been the recogni- 
tion that the morphologically normal majority of 
OB supergiants actually have some degree of N and 
He enhancements already, while the OBC minority 
have normal (i.e., main-sequence) abundances and 
the OBN group have more extreme mixing (Maeder 
& Conti 1994; Smith & Howarth 1994). The recent 
development of models for massive stellar evolution 
with rotation (Maeder & Meynet 2000) is an ad- 
vance toward realism since massive stars often rotate 
rapidly, but another dimension is added to the phase 
space that must be explored in order to compare ob- 
servations with theory. Fortunately, a new observa- 
tional diagnostic is added concurrently, namely the 
CNO abundances, since the degree of mixing also 
depends on the initial rotation. Indeed, progress 
has already been made toward understanding appar- 
ent discrepancies among the HRD locations of some 
mid-0 giants and supergiants in the Small Magel- 
lanic Cloud (SMC; Walborn et al. 2000), in terms of 
their CNO abundances and inferred initial rotations 
(Crowther et al. 2002; Hillier et al. 2003). Related 
results for main-sequence stars in the SMC giant H II 
region NGC 346 are reported by Bouret et al. (2003). 



Spectral classification of the hottest massive stars 
is relatively difficult, because of their scarcity and 
the paucity of reliable features in their optical spec- 
tra. Recently, an increased sample of high-quality 
data has supported classification developments for 
them, including introduction of the new types 02 
and 03.5 (Walborn et al. 2002). In this system, re- 
liance is placed on ratios of the selective emission 
lines (Walborn 2001) of N IV and N HI, primarily 
on an empirical basis. These new spectral types have 
yet to be modeled systematically and calibrated in 
terms of physical parameters, including confirmation 
that the 02 extension corresponds to higher effective 
temperatures as hypothesized. The analysis of these 
spectra is also difficult, because of uncertainties in 
the parameters and the atmospheric physics. E.g., 
the special behavior of the N IV and N III emission 
lines has to be reproduced and interpreted. Cur- 
rently, modeling and calibration developments are 
underway to address these issues; we shall present 
some relevant results here. 

Certainly additional spectral features to constrain 
the classification and analysis of the hottest spectra 
would be welcome. Drissen et al. (1995) pointed out 
certain little-observed N IV and O IV lines in the 
3400 A region of early- O spectra, that warranted 
further investigation in a larger sample including 
longer-wavelength classification standards. We have 
now accomplished such a survey, which will be fully 
reported by Morrell et al. (in preparation). Here 
we describe our unanticipated discovery of the first 
CNO anomalies in optical 02 spectra, which are rel- 
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evant to both the classification and modehng of these 
extreme objects. 



2. OBSERVATIONS 

The far-violet observations of the Large Magel- 
lanic Cloud (LMC) stars were carried out with the 
Cassegrain spectrograph at the Cerro Tololo Inter- 
American Observatory (CTIO) 1.5 m telescope dur- 
ing 2002 February 22-27. A 600 lines mm^^ grating 
blazed at 3375 A in second order combined with a 
Loral Ik CCD provided a 2 pixel resolution of 1.5 A 
and wavelength coverage from 3270 to 4180 A. The 
LMC stars discussed here have {V, B — V) rang- 
ing from (12.28, -0.23) for HDE 269810 to (13.68, 
-0.01) for LHlO-3061 (see Walborn et al. 2002 for 
detailed stellar parameters). Total exposure times 
ranging from 1 to 3 hours yielded signal-to-noise 
values between 60 and 150 per pixel at 3400 A, re- 
flecting the low (but relatively high at CTIO) at- 
mospheric transmission at that wavelength. The 
data for all but HDE 269810 have been smoothed by 
3 pixels for presentation here, and the radial veloci- 
ties have been corrected to rest. The data reductions 
and analysis were performed with standard IRAF^ 
routines running under Linux at the La Plata Ob- 
servatory. 

NGC '346-3 (13.50, -0.23) in the SMC was 
observed on the night of 2003 July 30-31 with 
the Boiler and Chivens spectrograph at the Las 
Campanas Observatory (LCO) Magellan I (Wal- 
ter Baadc) 6.5 m telescope. The detector was a 
2048 X 515 pixel Marconi CCD (pixel size 13.5 /um). 
A 600 lines mm~^ grating blazed at 5000 A pro- 
vided wavelength coverage from 3300 to 6400 A with 
a 2 pixel resolution of 3.1 A. The exposure time was 
10 min, producing a signal-to-noise of about 50 per 
pixel near both ends of the spectrogram but reaching 
300 at 5000 A. The data were reduced and analyzed 
with standard IRAF routines at LCO. 

Details of the blue-violet observations shown here 
can be found in Walborn et al. (2002). The ultravio- 
let observations discussed correlatively are from the 
Hubble Space Telescope (HST) Faint Object Spec- 
trograph (FOS) and Space Telescope Imaging Spec- 
trograph (STIS); they are fully described by Wal- 
born et al. (1995) and by Massey et al. (2004), re- 
spectively. 

As a check on the quantitative analysis of 
these modcratc-rcisolution data, an available high- 
resolution observation of HDE 269810 was also an- 
alyzed. This observation was made at the Euro- 
pean Southern Observatory (ESO) Very Large Tele- 
scope (VLT), with the UV- Visual EchcUe Spectro- 
graph (UVES) mounted at the UT2 8 m, on 2002 
November 29 under program 70.D-0164. The resolv- 

® IRAF is distributed by the National Optical Astronomy Ot 



ing power of this system is 70,000, providing a reso- 
lution of 0.05 A at the O IV A3400 lines. Two setups 
with several channels containing EEV and MIT/LL 
2k X 4k CCDs with 15 pixels covered wave- 
lengths from the far violet through Ha. Exposure 
times of 1200 s yielded S/N per 2 pixel resolution 
clement ranging from 80 at 3450 A to 150 at 4700 A 
and 120 at 6600 A. The data were reduced by the 
UVES pipeline under the Munich Image Data Anal- 
ysis System (MIDAS). 



3. RESULTS 

A sample of 27 Galactic and Magellanic Cloud 
early-0 stars has been observed in the 3400 A region; 
in most of them, the N IV and O IV features have 
comparable intensities and display systematic trends 
with spectral type and luminosity class, which will 
be discussed by Morrell et al. (in preparation). Sur- 
prisingly, however, the four LMC 02 III(f*) spectra 
in the sample separate into two groups, with O IV 
substantially stronger than N IV in one pair and 
the opposite in the other, while the SMC object is 
similar to the latter LMC pair, as shown in Fig- 
ure 1. This unexpected result prompted us to re- 
examine in a new light the blue- violet spectra of the 
same five stars shown by Walborn et al. (2002) and 
reproduced in Figure 2. It was immediately clear 
that the systematic differences between HDE 269810 
and Sk -68° 137 on the one hand, and LH64-16, 
LHlO-3061, and NGC 346-3 on the other, can be 
described consistently in terms of difFcrcnt relative 
line strengths of N vs. C. In particular, the first pair 
have weaker N III and N IV emission lines as well 
as N V absorption lines, along with the presence of 
C IV A4658, while all of the N lines are stronger in 
the latter group (except for the N V in NGC 346-3). 
Moreover, the He/H line ratios are seen to be larger 
in the latter group as well (although unclear in the 
SMC object due to nebular hydrogen emission and 
the lower resolution in the far violet). Presumably 
the O IV A4632, which is also prominent in the 
first pair, is blended with N HI A4634 in the other 
stars. The morphology strongly excludes ionization 
effects as an explanation and equally strongly sug- 
gests N/(C,0) and He/H abundance differences as 
the cause. 

Striking support for this conjecture is provided 
by the ultraviolet stellar-wind spectra of these 
stars. Those of HDE 269810, Sk -68° 137, and 
NGC 346-3 were discussed by Walborn et al. (1995) 
and are reproduced in Figure 3, along with that of 
LH64-16 obtained and analyzed by Massey et al. 
(2004). In HDE 269810, the C IV A1550 wind pro- 
file is comparable to that of N V A1240 and in 
Sk —68° 137 the C IV is actually somewhat stronger 
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than the N V, while in LH64— 16 N V is stronger and 

in NGC 346-3 much stronger than C IV. The UV 
spectrum of LHIO— 3061 has not yet been observed 
to our knowledge. 

Both the optical and UV spectra of LH64— 16 and 
NGC 346-3, as well as LHlO-3061 in the optical, 
relative to HDE 269810 and Sk -68° 137, are con- 
sistent with the characteristics of the ON class (Wal- 
born 1976, 1988, 2003). Hence, the first three spec- 
tra are hereby reclassified as 0N2 III(f*). This class 
has been interpreted in terms of mixing of CNO- 
cycled material into the atmospheres and winds of 
the stars (e.g., Maeder & Conti 1994). The presence 
of these effects at the earliest spectral types is a new 
development that may further challenge the models, 
however, as discussed below. 

Walborn et al. (2000) had already classified 
NGC 346-3 as 0N3 III(f*), purely on the basis of 
the very large N V/C IV wind-profile ratio in its 
UV spectrum. The ON designation was not re- 
tained when this star was moved into the 02 class 
by Walborn et al. (2002), partly because of concerns 
about a pure UV classification criterion and about 
the possibility that all such spectra might have simi- 
larly large N/C, but obvious only in the SMC object 
because the lower systemic metallicity allowed the 
very sensitive C IV feature to desaturate. However, 
LH64— 16 demonstrates that the latter concern is not 
valid, since the C IV can be weaker than N V at 
LMC metallicity. Thus, the proper classification of 
NGC 346-3 is also 0N2 III(f*), as confirmed here by 
the 3400 A observations. Haser et al. (1998) derived 
a very large N/C abundance ratio in NGC 346-3, but 
normal in Sk —68° 137, while Puis et al. (1996) gave 
normal He/H in both HDE 269810 and Sk -68° 137, 
consistent with the present classifications and inter- 
pretations. Bouret et al. (2003) have also derived 
an enhanced N abundance in NGC 346-3, as further 
discussed below. 



4. SPECTROSCOPIC MODELING 

In order to gain physical insight into the CNO line- 
strength differences outlined above, we have under- 
taken spectroscopic modeling of our blue and violet 
spectrograms. Example fits to the spectra of two 
stars are shown in Figure 4. 



4.1. Techniques 

Recent developments in model-atmosphere studies 
of hot stars routinely permit sophisticated metal- 
line blanketing in non-LTE to be considered, un- 
der either plane-parallel geometry (e.g., TLUSTY: 
Lanz & Hubeny 2003) or spherical geometry (e.g., 
CMFGEN: HiUier & Miller 1998). Bouret et al. 
(2003) have demonstrated the excellent consistency 



in synthetic spectra achieved between TLUSTY and 
CMFGEN in a combined UV and optical study of 
NGC 346-3, 0N2 III(f*). The former technique 
provides a more comprehensive treatment of photo- 
spheric metal-line blanketing, while the latter offers 
insights into wind properties. 

In the present application, we have utilized CM- 
FGEN, the current version of which is discussed by 
Hillier et al. (2003). The code solves the radiative 
transfer equation in the co-moving frame, under the 
additional constraint of statistical equilibrium. The 
model used here is similar to that applied by Bouret 
et al. (2003), except that some important additional 
ions are considered. In total, 27 ions of H, He, C, 
N, O, Si, P, S, Ar, Fe, and Ni are included, compris- 
ing a total of 2000 individual levels (grouped into 
800 superlevels), with a full array of 17,000 bound- 
bound transitions. The temperature structure is de- 
termined by radiative equilibrium. CMFGEN does 
not solve the momentum equation, so that a density 
or velocity structure is required. For the supersonic 
part, the velocity is parameterized with a classical 
/3-type law (with (3 = 1), which is connected to a hy- 
drostatic density structure at depth, such that the 
velocity and velocity gradient match at the inter- 
face. The subsonic velocity structure is set by the 
corresponding fully line-blanketed TLUSTY (v.200) 
model with log g = 4.0 (cgs). 

We have assumed a depth-independent Doppler 
profile for all lines when solving for the atmo- 
spheric structure in the co-moving frame, while in 
the final calculation of the emergent spectrum in 
the observer's frame, wc have adopted a radially- 
dependent turbulence, which reflects the effect of 
shocks due to wind instabilities. The minimum tur- 
bulence, ^, is set from the core strengths of H7 
and He II A4541. Incoherent electron scattering and 
Stark broadening for hydrogen and helium lines are 
adopted. We convolve our synthetic spectrum with 
a rotational broadening profile. Because of the in- 
termediate dispersion of our main observations, in- 
strumental effects prevent a reliable determination 
of wsini, so that 100 km has been adopted 
in each case, in accord with previous results for 
HDE 269810, Sk -68° 137, and NGC 346-3 (Puis 
et al. 1996). 

Temperatures of O stars are generally derived from 
spectroscopic fits to lines of He I-II (e.g., Conti 1973; 
Crowther et al. 2002). Since He I is absent or ex- 
tremely weak in 02 stars, we rely instead on the 
N IV- V ionization balance, i.e., N IV AA3479-85, 
4058 and N V AA4604-20. Of course, the nitro- 
gen abundance of an individual star is not known 
a priori, so that it is also a free parameter. With 
the stellar temperature set, we then derive the stel- 
lar radius from the absolute visual magnitude, for 
which we adopt values presented by Walborn et al. 
(2002). Since LHlO-3061 has not been observed in 



4 



the UV, we adopt a nominal wind terminal velocity 
of 3200 km s"-'^ for it, following prior UV analysis 
of Sk -68° 137 (Prinja & Crowther 1998); sources 
of terminal velocities for the other stars are given in 
notes to Table 1. Lacking the usual Ha wind diag- 
nostic for most of the stars, we estimate mass- loss 
rates from the strength of He II A4686, although we 
recognize that as a less reliable indicator. 

With regard to the helium abundance, gener- 
ally He/H = 0.1 (by number) is adopted imlcss 
noted otherwise. Potential nebular contamination 
of Balmer cores generally prevents accurate deter- 
minations of He/H from intermediate-dispersion ob- 
servations of O stars. For NGC 346-3, inspection 
of a high- resolution AAT /UCLES observation (Wal- 
born et al. 2000) supports He/H 0.1. In contrast, 
the CTIO 4 m observation of LH64-16 (Massey, 
Waterhouse, &: DeGioia-Eastwood 2000) analyzed 
here suggests a rather helium-rich chemistry, with 
He/H - 0.25. This same observation of LH64-16 
has been modeled independently by Massey et al. 
(2004), who find an even larger He/H abundance ra- 
tio. 

We estimate carbon and oxygen abundances from 
C IV A4658 and O IV AA3348-49, 3381-3410, 
respectively, while all remaining elements are set 
to 0.4 Zq {IMC) or 0.2 Zq (SMC), follow- 
ing, e.g., Russell & Dopita (1992). Solar abun- 
dances arc taken from Grevesse & Sauval (1998; 
log(C/H) + 12 = 8.52, log(N/H) + 12 = 7.92), ex- 
cept that the oxygen abundance has recently been 
revised downward to log(0/H) -|- 12 = 8.66 (Apslund 
2003). 



4.2. Results 

Early O giants in the Magellanic Clouds have pre- 
viously been studied quantitatively by Puis et al. 
(1996), who derived stellar temperatures of 55, 60. 
and 60 kK for NGC 346-3, Sk -68° 137, and 
HDE 269810, respectively, from plane-parallel pho- 
tospheric techniques, which neglected the effects of 
line blanketing. As with other recent redetermina- 
tions of 0-star temperature scales (e.g., Martins, 
Schaerer, & HiUier 2002; Crowther et al. 2002), we 
find 2.5-5 kK lower temperatures. Indeed, our re- 
sults closely mimic those of Bouret et al. (2003) for 
NGC 346—3, with stellar temperatures in the range 
52.5 < Tcff < 55 kK and luminosities in the range 
8 X 10^ < log(L/L0) < 2.2 X 10^ plus unclumped 
mass-loss rates of 1-2 x 10"^ Mg yr'^ (Table 1). 

More important in the present context arc elemen- 
tal abundances (Table 1). For the LMC 02 giants, 
we find only modest increases in the nitrogen abun- 
dances from the ISM ratios. In contrast, large N/C 
and N/0 ratios are obtained for the LMC 0N2 gi- 
ants, suggesting partially processed material at the 
stellar surfaces. For NGC 346-3 in the SMC, we 



confirm processed CNO abundances as derived by 
Bouret et al. (2003). Haser et al. (1998) previously 
estimated N > 0.4 Nq and C < 0.02 C© for this 
star from UV spectral fits. Relative to their corre- 
sponding ISM ratios, we find that N/C is enhanced 
in the atmospheres of LH64— 16, LHIO— 3061, and 
NGC 346-3 by factors of 46, 46, and 52, and N/0 
by factors of 14, 23, and 12, respectively. These fac- 
tors can be compared with those of order 100 derived 
by Crowther et al. (2002) and HiUier et al. (2003) for 
later-type O supergiants in the Magellanic Clouds. 



4.3. Uncertainties 

In order to estimate the uncertainties of our abun- 
dance results, we have undertaken additional mod- 
eling of two kinds. First, the fits to the data dis- 
cussed above have been re-examined, to determine 
the range of parameters, if any, that might yield 
comparable or perhaps even improved fits. For ex- 
ample, in the models upon which the results of Ta- 
ble 1 are based, the N IV A4058 emission line was 
well matched while the N V absorption lines were 
somewhat underpredicted. It was found that in- 
creasing the Teff of the two 0N2 stars, LH64— 16 
and LHIO— 3061, to 56 kK improves the fits to the 
N V lines but somewhat overpredicts the N IV emis- 
sion, as shown in Figure 4. Simultaneously, the fits 
to the very weak C IV A4658 emission line, the sole 
diagnostic of the carbon abundances in these data, 
were re-examined, with the result that C/Cq in the 
0N2 stars could be reduced to as low as 0.02. The 
only other possible changes indicated by this fur- 
ther modeling are small increases in O/O© for the 
two LMC 02 stars, HDE 269810 and Sk -68° 137, to 
0.6 and 0.7, respectively. In particular, no changes to 
the derived N abundances are indicated. To summa- 
rize, as best we can determine from the intermediate- 
dispersion data alone, the uncertainties in Teg are 
±1 kK, in the carbon abundances a factor of 2.5, and 
in the oxygen abundances a factor of 1.5. One could 
infer that the uncertainty in the nitrogen abun- 
dances is less than a factor of 1.5. 

A second, more powerful comparison is provided 
by the available VLT/UVES high-resolution spec- 
trogram of HDE 269810, with wavelength coverage 
extending into the yellow-red. Not only are the mul- 
tiple, weak O IV lines well resolved and defined in 
these data, additionally including AA3560— 63 and 
4632, but the strong C IV AA5801, 5812 emission 
lines are available to better constrain the carbon 
abundance, and N IV A6381 is also available. The 
good fit obtained with the same models used to 
analyze the intermediate-dispersion data is shown 
in Figure 5. The only significant discrepancies are 
in the O IV AA3404-12-14 triplet, which is pre- 
dicted in emission by models this hot but observed 
in absorption, and in higher members of the He II 
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Pickering series, which remain overpredicted. The 
improved fit to hydrogen and other heUum lines is 
noteworthy; the discrepancies at the lower disper- 
sions are believed to be caused primarily by line- 
broadening and resolution matching between the 
models and data. It is reassuring that the mass-loss 
rate based also on the Ha profile here agrees well 
with that determined previously from He H A4686 
alone. The stellar parameters that differ in the high- 
resolution analysis are Tofr, 55 kK; log L/_L0, 6.38; 
logM(M0 yr-i), -5.5; and Hc/H, 0.07; compari- 
son with the respective values in Table 1 indicates 
the reliability of that analysis. The derived C, N, O 
abundances relative to solar are 0.05, 0.6, and 0.8, 
respectively; i.e., the C abmidance is smaller by a 
factor of 4, while the N, O abundances are larger 
by respective factors of 1.5 and 2. It follows from 
these results that HDE 269810 already has a an ap- 
preciable N/C enlianccmcnt over the LMC ISM, by 
a factor of 23. 

The uncertainties in the solar CNO abundances 
given by Grevesse & Sauval (1998) are ±0.06 dex, 
or a factor of 1.15. The corresponding Magellanic 
Cloud ISM abundance uncertainties from Russell & 
Dopita (1992) are 0.13-0.18 dex for C, 0.15-0.20 dex 
for N, and 0.06-0.10 dex for O, or average factors 
of 1.4, 1.5, and 1.2, respectively. Thus, the uncer- 
tainties in the ratios are dominated by the stellar 
abundances, which we have estimated above to be 
factors of 2.5 4 for C, and 1.5-2 for N and O. The 
worst case is the N/C ratio, which in principle could 
be uncertain by a factor of 6. However, such a factor 
cannot apply randomly to the results, since it could 
invert the differential spectral morphology, which is 
not reasonable; on the other hand, a systematic ef- 
fect of that magnitude cannot be ruled out with the 
present data. It is interesting that all of the C rede- 
terminations are lower, while the N and O are some- 
what higher. Thus, the N/C ratios may be even 
larger than given in Table 1. Clearly it is highly 
desirable to reanalyze all of these stars with high- 
resolution data, at the same time determining their 
mass-loss rates from Ha and UV data, to refine the 
quantitative results. 

5. DISCUSSION: EVOLUTIONARY 
INTERPRETATIONS 

Curiously, the luminosities, and hence the masses, 
estimated by Walborn et al. (2002) for the 02 III(f*) 
stars HDE 269810 (150 Mg, average of two meth- 
ods) and Sk -68° 137 (100 Mq) are considerably 
higher than those for the 0N2 III(f*) stars LH64-16 
(62 Mq), LHlO-3061 (81 AIq), and NGC 346-3 
(82 Mq). The opposite relationship might have 
been expected, since both mass-loss rates and mix- 
ing increase with mass (Meynet & Maeder 2000). 
On the other hand, perhaps the longer evolutionary 



timescales at the lower masses favor the observabil- 
ity of mixing effects. 

Until recently, these abundance results would have 
been interpreted simply in terms of more advanced 
evolutionary ages for the 0N2 than the 02 objects, 
which remains a possibility. (Of course, mass trans- 
fer in a binary system is an alternative mechanism 
likely responsible for at least some members of the 
ON class — see Bolton & Rogers 1978; Levato et al. 
1988; Walborn & Howarth 2000; Lennon 2003). In 
the context of evolutionary models with rotation for 
massive stars (Maeder & Meynet 2000; Meynet & 
Maeder 2000), however, the degree of mixing of pro- 
cessed material to the surface at a given evolutionary 
age is a function of the initial rotational velocity. In- 
deed, CNO abundances in massive stars now become 
critical diagnostics of their initial parameters and 
evolutionary histories. For a given metallicity, the 
observed Toh , L, and CNO abundances in principle 
enable derivation of the corresponding initial mass, 
current age, and initial rotational velocity. Com- 
plete grids of evolutionary tracks for the most mas- 
sive stars, at Magellanic Cloud metallicities and with 
a range of rotational velocities, are currently under 
development in Geneva. Meanwhile, we display il- 
lustrative comparisons in Figures 6-8, courtesy of 
G. Meynet and A. Maeder. 

Figure 6 is a theoretical HR Diagram in which 
the locations of the stars discussed here are com- 
pared with Geneva evolutionary tracks for initial 
masses of 60 and 120 Mq, at both LMC and SMC 
metallicities (mass fractions Z = 0.008 and 0.004, 
respectively), and with an initial equatorial rota- 
tional velocity of 300 km s^^. All of the stars lie 
near the zero-age main sequence (ZAMS), which is 
consistent with an early evolutionary stage for the 
two objects with normal spectra; the present spec- 
troscopic analysis again indicates higher masses for 
them than for the 0N2 objects. However, such a 
location is surprising for the latter, highly mixed 
stars. Their initial masses still place them above the 
Humphreys-Davidson Limit, so they cannot be in- 
terpreted as post-red supergiants evolving blueward. 
A more relevant possibility may be that of "homo- 
geneous" evolution, along tracks that rise and curve 
blueward from the ZAMS in the HRD, which is pre- 
dicted in the cases of strong mixing and/or very high 
initial rotational velocities (Langer & El Eid 1986; 
Maeder 1987; Langer 1992; Meynet & Maeder 2000)! 
The 60 Mq homogeneous track at solar metallic- 
ity (Z ~ 0.020) with an initial rotational velocity 
of 400 km s^^ from the last reference is also repro- 
duced in Fig. 6; it trends near the locations of the 
0N2 III(f*) points. Indeed, Bouret et al. (2003) 
have concluded that such a track is required to un- 
derstand NGC 346-3 as coeval with the other early 
O stars in the cluster. 
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Figures 7 and 8 compare our derived N/C and N/ O 
ratios, respectively, relative to their corresponding 
interstellar media, taken at face value from Table 1 , 
with the predictions of the same Geneva models. 
The N/C ratios in the 0N2 III(f*) stars are the 
largest reported to date for near-main-sequence ob- 
jcc;ts and are substantially larger than the standard 
model predictions, suggesting even more rapid mix- 
ing and perhaps providing additional support for ho- 
mogeneous evolution. In contrast, the N/0 ratios 
are compatible with the standard model predictions 
and indicate early operation of the ON cycle, as ex- 
pected in only the most massive stars. 

On the other hand, the even larger N enhance- 
ments derived for later-type O supergiants in the 
Magellanic Clouds by Crowther et al. (2002) and 
Hillier et al. (2003) may imply that the 0N2 III(f*) 
stars are simply predecessors of those objects along 
redward evolutionary tracks. Again, much stronger 



mixing effects near the main sequence than predicted 
by current models, at least for N/C, are then re- 
quired by our observations. Further analysis of these 
phenomena, with improved data and relative to the 
latest models for rotating massive stars, will clarify 
these issues and provide an improved understanding 
of their post-main sequence evolution. 

We are very grateful to Georges Meynet and Andre 
Maeder for their substantial contributions to the dis- 
cussion of these results, including Figures 6 8. We 
also thank the STScI Director's Discretionary Re- 
search Fund for publication support. PM's partici- 
pation was supported by NASA through grant GO- 
8633 from the Space Telescope Science Institute, 
which is operated by AURA, Inc., under contract 
NAS5-26555. We thank an anonymous referee for 
useful recommendations. 
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(a) Walborn et al. 1995; (b) Prinja & Crowther 1998; (c) Massey et al. 2004; 
(d) Bouret et al. 2003; (e) Walborn et al. 2002 
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3300 3400 3500 3600 3700 3800 3900 4000 4100 



Wavelength (A) 

Fig. 1. — Rectified linear-intensity, far-violet spectrograms of five 02 giants in the Magellanic Clouds, 
separated by 0.5 continuum units. The absorption lines identified below are O IV AA3348— 49, 3381—85, 

3397, 3410-14; N IV AA3479-85 (the latter a blend of two lines); interstellar Ca II A3933; He A3970; He II 
A4026; and H(5 A4101. The N IV A4058 emission line is identified above (and the O IV, N IV absorption- 
line brackets are replicated). In the spectrogram of LHIO— 3061, an ovcrsubtractcd [O II] A3727 nebular 
emission doublet has been patched out. Note the different relative strengths of the O IV and N IV features 
between the 0N2 and 02 classes. 
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Fig. 2. — Blue-violet spectrograms of the same stars and with format as in Fig. 1. The absorption 

lines identified below are He A3970, RS A4101, and H7 A4340; He II AA4026, 4200, 4541, 4686; and N V 
AA4604-20. At the top, N IV A4058; Si IV AA4089-4116; and N HI AA4634, 4640-42 emission lines are 
identified. Above HDE 269810, O IV A4632, Si IV A4654, and C IV A4658 emission lines are identified; 
this O IV line had been identified with Si IV A4631 by Walborn et al. (2002), but we have now definitively 
established from its wavelength and intensity in high-resolution data (Fig. 5) that the O IV contribution 
is dominant. Hydrogen nebular emission lines have been truncated in NGC 346—3. Note the different 
N/C and He/H line ratios between the 0N2 and 02 classes (the latter as permitted by nebular H emission 
contamination) . 
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N V IV V Si IV C IV He II N IV 



1200 1300 1400 1500 1600 1700 

Wavelength (A) 

Fig. 3. — HST far-ultraviolet spectrograms of four stars shown in the previous figures. The longer ordinate tick marks give the zero levels and 
the continua are normalized to 1.0. The spectral features identified below are N V AA1239-43, O IV AA1338-43, O V A1371, Si IV AA1394-1403 

(interstellar), C IV AA1548-51, He II A1640, and N IV A1718. Note the difference between the relative strengths of the N V and C IV features 
between the upper and lower pairs of spectra; the O IV and O V features also appear weaker in the 0N2 spectra. 
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Fig. 4. — Atmospheric model fits {dotted lines) to key diagnostic features in two of the intermediate- 
resolution spectrograms {solid lines). The 56kK fit to LH64— 16 is shown (see text). The spectral lines iden- 
tified in the top panel arc O IV AA3348-49, 3381-85-90-97-3410, 3404-12-14 and N IV AA3479-83-85. 
In the middle panel, N IV A4058 is identified, and in the bottom panel, N V AA4604-20, O IV A4632, and 
C IV A4658. 
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Fig. 5. — Atmospheric model fit {dotted lines) to the VLT/UVES high-resolution observation of 
HDE 269810 {solid lines). The spectral lines identified by panel from top to bottom are O IV AA3348— 49, 
3381-85-90-97-3410, 3404-12-14, 3560-63 and N IV AA3479-83-85; N IV A4058; N V AA4604-20, 
O IV A4632, and C IV A4658; C IV AA5801-12; and N IV A6381. 
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Fig. 6. — Theoretical HR Diagram comparing the observed stars with Geneva evohitionary tracks at both 
LMC and SMC metalhcities, for an initial rotational velocity of 300 km s^^. Also included is a 60 Mq 
homogeneous track at solar metallicity, for an initial rotational velocity of 400 km s~^. Squares, LMC stars; 
triangle, SMC star; solid lines, LMC tracks; short dashed lines, SMC tracks; long dashed line, homogeneous 
track. 
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Fig. 7. — Derived stellar N/C abundance ratios relative to their initial (interstellar) values compared with 
predictions from the models in Fig. 6. The 60 Mq tracks are toward the left and the 120 Mq toward the 
right; short vertical line segments on the former indicate the end of the hydrogen-burning main sequence. 
Other symbols and line styles as in Fig. 6. 
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Fig. 8. — Derived stellar N/0 abundance ratios relative to their initial (interstellar) values compared with 
predictions from the models in Fig. 6. Format, symbols, and line styles as in Fig. 7. 
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